Comment on "Can Disorder Induce a Finite 
Thermal Conductivity in ID Lattices?" 

In a recent letter, Li et al [Q have reported that the 
steady state of the disordered harmonic chain is not 
unique and depends on initial conditions. Their claim 
is based on a molecular dynamics simulation using Nose- 
Hoover thermostats to model the heat baths. We point 
out that the uniqueness of the nonequilibrium steady 
state of the disordered harmonic chain for a large class 
of heat baths has been proven exactly . In this note 
we consider a particular case where it is easy to explicitly 
demonstrate the uniqueness of the steady state. This is 
the case where the heat bath is modelled by Langevin 
dynamics. The results of Li et al appear to be either due 
to insufficient equilibration times or an artifact of using 
Nose-Hoover thermostats. 

Consider heat conduction through a one-dimensional 
disordered harmonic chain. Particles i = 1,2...N with 
random masses are connected by harmonic springs of 
equal strengths. The Hamiltonian of the system is thus 
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where {xi} are the displacements of the particles about 
their equilibrium positions, {pi} their momenta and {mi} 
are the random masses. We put the boundary conditions 
xq = xn+1 = 0. The particles in the bulk evolve through 
the classical equations of motion while the boundary par- 
ticles, namely particles 1 and N are coupled to langevin 
heat baths as in Ref. 

Thus we consider the following equations of motion for 
the particles: 

mix'i = -2xi + X2- xi + r]L{t) 
mxi = -{2xi - xi^i - xi+i) I = 2, 3...{N - 1) 
ruNXN = -2xN + XN-i) -±N +riRit), (2) 

where 77^ and 77/? are gaussian white noises with the cor- 
relations {'qL,R{t)-qL,R{t')) = 2TL,R5{t-t'). Denoting the 
coordinates and momenta collectively by the variables qi 
so that {gi,g2---92iv} = {xi,X2..-xn ,Pi,P2---Pn} Eqs. (||) 
can be written in the form: 
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where the vector 77 has all elements zero except jyAr+i = 
rjL and •q2N = "HR and the 2N x 2N matrix a is given by: 


with 


J 

Jki — ~^k,i/''nk; ^ki = 25fc^, — — Skd+i 

Iki = 5k,i{5k,i/mi + Sk,N/mN). 
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In the steady state {d{qkqi) / dt) = 0. From this and using 
Eq. (H) we get the matrix equation 


where b is the correlation matrix with elements bki = 
{qkqi) and dki = dkji2TLSk,N+i + 2T7j4,2Ar)- We can 
invert this equation to obtain b and thus all the mo- 
ments including the local temperatures Ti = pf /mi. The 
uniqueness of the steady state then depends on whether 
or not Eq. (|) has a unique inverse. For chains of finite 
length and for given disorder realizations it is easy to ver- 
ify numerically that Eq. does have a unique inverse. 
We also find that a molecular dynamics simulation (using 
a simple Euler discretization), with langevin heat baths, 
reproduces the exact temperature profile and is indepen- 
dent of initial conditions. This is shown in Fig. (|l|). For 
the iV = 20 lattice, averaging over 10^ time units is suffi- 
cient to achieve steady-state values. The = 40 data is 
averaged over 10* time units. We find that equilibration 
times increase rapidly with system size. 

Finally we note that for Langevin heat bath dynam- 
ics, the unique steady state distribution is infact known 
exactly H and is a gaussian given by: 
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FIG. 1. The exact temperature profiles as obtained from 
inverting Eq. Q are compared with those from molecular 
dynamics simulations for two lattice sizes A'^ = 20 and A'' = 40 
(inset). 
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